Nonequilibrium Kondo Effect in a Multi-level Quantum Dot near singlet-triplet 

transition 
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The linear and nonlinear transport through a multi-level lateral quantum dot connected to two leads 
is investigated using a generalized finite-(7 slave-boson mean field approach. For a two-level quantum 
dot, our calculation demonstrates a substantial conductance enhancement near the degeneracy point 
of the spin singlet and triplet states, a non-monotonic temperature-dependence of conductance and 
a sharp dip and nonzero bias maximum of the differential conductance. These agree well with 
recent experiment observations. This two-stage Kondo effect in an out-of-equilibrium situation is 
attributed to the interference between the two energy levels. 

PACS numbers: 72.f5.Qm, 73.23.Hk, 73.40.Gk, 73.50.Fq 
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Since the observation of the Kondo effect in semicon- 
ductor quantum dot (QD) with odd electron number [1], 
there have been a great deal of experimental and the- 
oretical investigation into this many-body phenomenon. 
With easy control of major parameters of these artificial 
atoms in a wide range, QDs facilitate exploration of the 
Kondo physics at very different environments. As a con- 
sequence, new types of Kondo effects have been discov- 
ered in experiments. Very recently, a surprising Kondo- 
enhanced conductance has also been detected for an even 
number of electrons in both vertical [2] and lateral [3,4] 
QD configurations, which is markedly different from the 
conventional spin- 1/2 Kondo effect. 

A real QD contains more than one energy levels and 
usually two electrons incline to occupy the lower level, 
resulting in a singlet local spin state. Decreasing the 
spacing of the two energy levels by applying a magnetic 
field makes it possible for two electrons to occupy two 
different levels due to the exchange interaction J, form- 
ing the local spin-triplet state. Scaling theory [5,6] and 
numerical renormalization group analysis [7,8] revealed 
that this singlet-triplet transition gives rise to an anoma- 
lous enhancement of conductance at low temperature. 

The physics in the lateral and vertical QD configu- 
rations are different. Recent observation on the lat- 
eral QD has disclosed a two-stage Kondo effect in the 
temperature-dependent conductance and differential con- 
ductance, and favored an interpretation in terms of a sin- 
gle conduction channel in two leads. [4] So far, there have 
been only a scaling analysis [6] and a numerical renormal- 
ization group calculation [8] on this problem in the low 
bias regime. Its clear physical picture and the nonlin- 
ear conductance at finite bias need further exploring. In 
this letter we investigate the singlet-triplet Kondo effect 
in a lateral QD having two energy levels in an out-of- 
equilibrium situation, by generalizing the finite-?/ slave- 
boson mean-field (SBMF) approach, which was initially 
developed by Kotliar-Ruckenstein [9] and extended by 
us to study Kondo-type transport through QD [10] and 
coupled QDs [11]. Our results manifest not only a sub- 



stantial enhancement of conductance near the local spin 
singlet-triplet transition, but also a two-stage Kondo ef- 
fect in temperature-dependent conductance and differen- 
tial conductance. 

The Hamiltonian of the QD with two orbital energy 
levels (j = 1,2) connected to two leads (77 = L/R) 
can be written as H = Hl + Hr + Ho + Ht 1 in which 



H.c.) describe the left and right leads and the tunneling 
between each lead and the QD (assuming a single con- 
duction channel available per lead), respectively, and Hrj 
represents the isolated QD (a, /3, 7,er = ±1): 
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where Cj in the first term is the single particle energy, 
the second and third terms denote the intra- and inter- 
level Coulomb interactions, and the last term gives a 
ferromagnetic exchange coupling J > due to Hund's 
rule. The energy spacing 8 = e 2 — £i can be controlled 
by an external magnetic field. However, considering the 
small <?-factor in, for example, GaAs, we neglect the Zee- 
man splitting in the QD, as done in previous theoretical 
treatment. [7,8] To keep the discussion simple we assume 
Uj = U 12 = U and Vjr, = V {j = 1 , 2 and 77 = L , R). 

Following the SBMF scheme, [9] we introduce six- 
teen auxiliary Bose fields associated with sixteen single- 
particle eigenstates of the isolated multi-level QD Hamil- 
tonian (1), as shown in Table I. For example, boson oper- 
ator d\s z is with the spin-triplet state having spin (1, S z ), 
and dgi is with the spin-singlet state | fj,, 0) (it is the low- 
est energy state in all three spin-singlet states) . Note that 
the splitting of energies between this spin-singlet state 
and spin-triplet state A = J — S can be tuned by a mag- 
netic field. Therefore, for 6 — J the three triplet states 
diSz an d the singlet state doi are degenerate. Necessarily, 
the completeness relation for these slave-boson operators, 
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the correspondence between fermions and bosons, n JCT 
t „ — n fn. = Ji n _i_^t 
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be imposed to confine the enlarged Hilbert space. More- 
over, in the combined fermion-boson representation, the 
QD fermion operators cj CT and Cj a in the hopping term 
are exppressed as z\ _ct and c 1CT z, CT , respectively, where 

jo jo J J 

Zj a consists of all the boson operator sets which are as- 
sociated with the physical process that a a-spin electron 
of the jth level is annihilated (j + 1 j, a ^ a): 

z ja = QJ a 1/2 (1 - Q^)- 1 ' 2 \jp ja +d a d la 
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TABLE I. Sixteen Eigenstates, spin quantum numbers 5, 
S z and energies E for the isolated QD with two levels 1 and 
2, and the assigned slave-boson (SB) operators. 
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Within the slave-boson scheme, the Hamiltonian of the 
system can be replaced by the following effective Hamil- 
tonian in terms of auxiliary boson operators plus the con- 
straints incorporated via the Lagrange multipliers A and 
\j a : 

rf,k : (T j,(7 
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From the effective Hamiltonian (3) one can derive 
equations of motion of the slave-boson operators. Then 
we use the mean-field approximation in these equations 
and in the constraints, in which all the boson operators 
are replaced by their expectation values. The negligible 
Zeeman splitting helps to reduce the number of the inde- 
pendent expectation values, such that we have pj a = pj, 
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and Xj a = Xj. With 



the help of the Langreth analytical continuation rules 
for the close time-path Green's function (GF), [12] these 
equations of motion can be closed in terms of the dis- 
tribution GF Gf(uj) of the QD, leading to the following 
self-consistent set of equations = 1, 2): 
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'-Vj + 2(A - 2Ai - 2A 2 + 4(7 - 2J)f = 0, (12) 



in which V j = ^ l J dujGf (w) (w - ij) and 
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with 



V{uj) = (w - ei)(u - e 2 ) ± iTi(u - ?i) ± iT 2 {uJ - e 2 ), 

(14) 
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e„fe) being the coupling constant between the QD and the 
lead. f n (u>) is the Fermi distribution function of 77th lead 
(77 = L/R). The current J through the two- level QD can 
be divided into three parts: the contributions from the 
ground level I\, from the second level I 2 , and the current 
Ii resulting from the interference between the two energy 
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levels. Straightforwardly, we have J = I\ + J 2 + h = 
X J-™ [Ti M + W + [f L (w) - , in which 

Ti(w) = fi(w - e 2 ) 2 |X>(u;)|~ 2 , (15) 
T 2 (w) = fa(w - ei) 2 |D(w)|~ 2 , (16) 
Ti(w) = 2f if 2 (w - e 2 )(w - ei)|D(w)|- 2 (17) 

are the transmission probabilities. 

In the linear limit the conductance can be written as 

G = ^[T l( 0) + T 2 (0) + 11(0)] = ^ {flh + f ae " l)S 



P(0)| 2 
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In numerical analysis, we fix the parameters of the QD: 
U = 4, J = 2 with ei = —2.5 or —3, to guarantee nearly 
two electrons dwelling in the QD (for — 4 < A < 1.2, 
1.45 < N < 1.6 at ei = -2.5; 1.5 < N < 1.7 at ei = -3), 
then focus our attention on the A-dependent Kondo ef- 
fect in transport. Here T is always taken as the energy 
unit and ej is measured from the Fermi level of leads. 
Fig. 1 and the inset figure show the total linear conduc- 
tance G and its three parts as functions of the energy 
splitting A between the spin-singlet and -triplet states in 
the QD at zero temperature. There appears a peak in 
G near the singlet-triplet degeneracy point A = 0. In 
the spin-singlet-state-dominated regime (singlet regime) , 
A < 0, the contribution of the second level, G 2 , and that 
of the interference between the two levels, Gi, are almost 
zero. Only the first level is the active conduction chan- 
nel. With increasing A (decreasing S), the second level 
starts to carry current, such that G rises up. Around the 
degeneracy point, the interference effect appears and its 
contribution G; is negative. The rapid enhancement of 
the interference effect causes an abrupt decline of G in 
the spin-triplet-state-dominated regime (triplet regime) 
A > and a sharp peak in G near the transition point. 
Note that the nonzero value of G in the singlet regime 
is determined by the position in the Coulomb blockade 
valley due to N < 2 in our model QD parameters. 

Fig. 2 illustrates the behavior of the conductance G as a 
function of temperature up to 0.3. It is believed that the 
SBMF approach is correct for describing the Kondo ef- 
fect in strongly correlated systems at temperatures lower 
than the Kondo temperature Tk- For the systems un- 
der consideration, the Kondo temperatures Tk are esti- 
mated to be about 0.21 and 0.32 if only the first level 
is considered. [13] Therefore we can cautiously believe 
that the present calculation gives a rational description 
for the temperature dependence of the Kondo contribu- 
tion to the conductance. In the cases A = and 0.5, the 
conductance rises monotonously with decreasing temper- 
ature. In the cases A = 1 and 1.2, however, a "hump" 
behavior appears. After riseup with lowering tempera- 
ture at first, G reaches a peak at a certain value of T be- 
fore decreases, indicating a two-stage Kondo effect in this 



triplet regime. These results arc in qualitatively agree- 
ment with the scaling analysis of Pustilnik and Glazman 
[6] , where a two-stage Kondo effect for an S > 1 ground 
state (i.e., the triplet regime in this letter) was suggested. 
Moreover, we can see from Fig. 2 that more pronounced 
hump appears in deeper triplet regime (bigger A). It is 
also worth mentioning that, for moderate A, for example, 
A = 1 at ei = -2.5 [Fig. 2(a)] and A = 0.5 at e 1 = -3 
[Fig. 2(b)], the conductance G may slightly increase at 
very low temperatures, exhibiting a "shoulder" behavior, 
which satisfies very well with the recent measurements in 
Ref. [4] (inset of Fig. 2). 

To simplify the calculation of nonlinear transport un- 
der a finite bias voltage between the two leads, we assume 
a symmetric voltage drop, = —fiR = eV/2. The cal- 
culated zero-temperature current J as a function of bias 
voltage V in the case of U = 4, J = 2 and ei = —2.5 
are shown in Fig. 3(a) for A > 0.5. With increasing the 
bias voltage in the range shown in the figure, 1) the cur- 
rent / increases monotonically for A < 0.5; 2) / ap- 
proaches a peak value followed by a drop for A > 0.5; 
3) the peak current appears at smaller bias for bigger 
A. Correspondingly the voltage-dependent differential 
conductance dl/dV [Figs. 3(b)] shows a non-zero-bias- 
maximum at triplet regime A > 0.5. While with de- 
creasing A, this sharp dip structure gradually transits to 
a conventional out-of-equilibrium Kondo behavior, i.e., 
zero-bias-maximum, at A = 0.5. As shown in the in- 
sets in Figs. 2 and 3(b), where the separate contributions 
of the first and second levels as well as the interference 
term, to G and dl/dV are shown respectively, the hump 
behavior and the shoulder structure in G-T curves, and 
the sharp dip in dl/dV-V curves, can all be attributed 
to the strongly weakening of the interference effect due 
to increasing temperature or bias voltage. 

In conclusion, we have generalized the finitc-C/ SBMF 
approach to the case of multi- level QD, and employed 
it to investigate the singlet-triplet Kondo phenomena 
in linear and nonlinear transport through a two levels 
QD in the N ~ 2 regime. Our investigation revealed a 
two-stage Kondo effect, a non-monotonic temperature- 
dependent conductance and a sharp dip structure in dif- 
ferential conductance, depending on the energy splitting 
A between the spin-singlet and triplet states. These are 
in agreement with recent experimental observation. The 
predicted behavior of linear conductance is in consistent 
with existing theoretical analysis. An interpretation for 
this two-stage Kondo behavior in G-T and dl/dV-V has 
been provided in terms of the interference between the 
two levels. 
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Figure Captions 

Fig.l Linear conductance at the absolute zero tempera- 
ture as a function of the splitting between the spin-singlet 
and triplet states A in QD of U = 4 and J = 2 with 
ei = —2.5 and —3. Inset: Corresponding contributions 
of three parts Gi, G2 and Gi for ei = —3. 

Fig. 2 Temperature dependence of linear conductance 
with A = 1.2, 1, 0.5, and for (a) ei = -2.5 and (b) 
ei = —3. Insets: Corresponding contributions of three 
parts Gi, G 2 and G; for (a) A = 1 and for (b) A = 1.2. 

Fig. 3 (a) Zero-temperature I-V characteristic curves for 
QDs with U = 4, J = 2, and e x = -2.5 at A = 1.2 - 0.5. 
(b) Zero-temperature differential conductance dl/dV as 
a function of the external voltage at A = 0.5, 1.0, and 



1.2. Inset in (b): Corresponding contributions of three 
parts dh/dV, d! 2 /dV and dl-JdV. 
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